Acoustic plasmons and "soundarons" in graphene on a metal gate 
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We demonstrate that single-layer graphene in the presence of a metal gate displays a gapless 
collective (plasmon) mode that has a linear dispersion at long wavelengths. We calculate exactly 
the acoustic-plasmon group velocity at the level of the random phase approximation and carry out 
microscopic calculations of the one-body spectral function of such system. Despite screening exerted 
by the metal, we find that graphene's quasiparticle spectrum displays a very rich structure charac- 
terized by composite hole-acoustic plasmon satellite bands (that we term for brevity "soundarons" ) , 
which can be observed by e.g. angle-resolved photoemission spectroscopy. 



Introduction. — Graphene, a monolayer of carbon 
atoms packed in a two-dimensional (2D) honeycomb lat- 
tice, is a gapless semiconductor whose recent isolation has 
generated an enormous deal of interest in various fields 
of condensed matter physics pQ. Electrons in graphene 
move as if they were massless Dirac fermions (MDFs) 
with a velocity v which is roughly three-hundred times 
smaller than the velocity of light in vacuum. Apart from 
a few notable exceptions [2], until 2009 many electronic 
and optical properties of graphene could be explained 
within a single-particle picture in which electron-electron 
(e-e) interactions are completely neglected. The discov- 
ery of the fractional quantum Hall effect in graphene [3] 
represents an important hallmark in this context. By now 
there is a large body of experimental work [4HS] showing 
the relevance of e-e interactions in a number of key prop- 
erties of graphene samples of sufficiently high quality. 

Recently, it has been demonstrated that the MDF 
energy spectrum of graphene is significantly altered 
by many-particle effects: robust "plasmaron" satel- 
lite bands have been first predicted [TJ [8] and then 
observed [5 in quasi-freestanding graphene grown on 
hydrogen-terminated SiC. Plasmarons [9 are composite 
quasiparticles that emerge in an electron liquid when the 
interactions between electrons and plasmons - the collec- 
tive density oscillations of the electron liquid [10 - are 
particularly strong. A deep understanding of the cou- 
pling between charge carriers and plasmons in graphene 
is highly desirable because of the large potential of this 
material [11 in the context of plasmonics [T2] . 

This Letter is motivated by the large amount of work 
that has been carried out by many groups around the 
world [13] to grow high quality graphene sheets on 
transition-metal templates such as Ru, Ir, Ni, Pt, Cu, 
and Au. Depending on the metal, there can be weak 
or strong hybridization between graphene 7r and metal d 
bands. In the former case, electrons close to the Fermi 
energy in the graphene sheet are still described by a MDF 
Hamiltonian: the main qualitative role of the metal is to 
screen e-e interactions between MDFs, weakening them. 
In this Letter we present a careful analytical and numeri- 
cal analysis of the impact of a metal gate on the collective 



behavior of the electron gas in a nearby graphene sheet: 
we show that, contrary to expectation, screened e-e in- 
teractions in graphene on a metal gate (G/M) can lead to 
intriguing many-body effects and rich spectral features. 

Model Hamiltonian and acoustic plasmons. — We 
model the electron gas in G/M by the following contin- 
uum Hamiltonian (hereby written for a single channel): 

H = hv V^L(^/3 • k)$ h j + — v d(v)pqP-q • 

k,a,P q^O 

(i) 

Here v is the MDF velocity, S is the sample area, and 
Pq = Y,k,a ^k-q,a^k,a is the density operator. Greek 
letters are honeycomb-sublattice-pseudospin labels and 
a = (o- x ,cr y ) is a 2D vector of Pauli matrices. Since 
we are interested only in the role of screening played by 
the metal gate, in writing Eq. we have neglected hy- 
bridization i) between the graphene sheet and the metal 
and ii) between the surface plasmons of the metal [T4] 
and the low-energy "sheet" plasmon of the graphene elec- 
tron gas. The metal is thereby modeled as a grounded 
conductor parallel to the graphene sheet. The function 
Vd(q) represents the Fourier components of the e-e inter- 
action, altered by the presence of the grounded gate [15] : 
Vd(q) — 27re 2 [l — exp(— 2qd)]/q, d being the graphene- 
metal gate distance. Due to screening exerted by the 
metal, Vd(q) is regular at q = for any finite d, i.e. 
Vd(q —> 0) = 47re 2 d. In the limit d — » oo one recovers the 
usual 2D Fourier transform of the long-range Coulomb 
potential, i.e. Voo(q) = 27re 2 jq. 

Whether or not the model represented by Eq. ([TJ is ad- 
equate to describe the electron gas in epitaxial graphene 
grown on metal substrates [13] (the metal being one of 
those that hybridizes weakly with graphene such as Pt 
or Cu) can only be determined a posteriori by check- 
ing the predictions of this Letter against experiments. 
Our results are also relevant, though, for mechanically- 
exfoliated flakes in the proximity of a metal gate. 

The collective density excitations of the Hamiltonian T-L 
are completely determined by the density-density linear- 
response function, Xpp(q^) = {{f>q\ P-q))w I f S, with 
((A^B))^ the usual Kubo product [ID] . Within the ran- 
dom phase approximation (RPA) [ID] [16] XppiQi 00 ) ls 
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where Xo(q^) 1S the well-known [17 response function 
of a 2D noninteracting gas of MDFs at arbitrary doping 
n. For future purposes, we also introduce the coupling 
constant a ee = e 2 /(hv), the Fermi wave number kp = 
yj4n\n\/(g s g w ), g s = 2 (g v = 2) being a spin (valley) 
degeneracy factor, and the Fermi energy sp = vkp. The 
collective modes of the Hamiltonian H are particle-hole 
symmetric since Xo(q^) is so. 

The dispersion uj ac = oj ac (q) of the collective (plas- 
mon) mode can be found by solving the complex equa- 
tion e(q,(jj) = in the regions of the q,uj plane where 
^sm [xo(q^)] — 0- Since Vd(q) is regular at q = 0, we 
expect a gapless acoustic plasmon, u ac (q — 0) = c s g, 
rather than the usual ("unscreened") plasmon [7J [17] 
oc yfq for q — >• 0. The latter behavior is formally re- 
covered only at d = oo. We now proceed to derive an 
exact expression for the RPA group velocity c s . Fol- 
lowing Santoro and Giuliani [18], we first introduce the 
power expansion LO ac (q) = c s q + C2q 2 + C3# 3 + . . . , for the 
acoustic-plasmon dispersion, and then define the function 
F(q) = e(q, c s q+C2q 2 -\-c3q 3 + . . . ). In the long- wavelength 
limit -F(#) admits the following Laurent- Taylor expansion 
F{q -^0) = /o + /ig + /2^ 2 + -" , where the coefficients 
fi can be derived from the analytical expression [17] for 
Xo(q,u). For s(<7,cj) = to be valid we have to require 
that the coefficients fi vanish identically. The coefficient 
fo depends only on c s and by equating its expression to 
zero we arrive after some tedious but straightforward al- 
gebra at the desired expression for the acoustic-plasmon 
group velocity c s : 



/ A(a ee dk-F) 
A(a ee dk^) — 1 



(3) 



with A (x) = [l+l/(2g s g v x)} 2 > 1. Note that the "sound" 
velocity c s is larger than v for any value of the dimension- 
less product a ee dkp. This implies that, within RPA, an 
undamped acoustic plasmon is always present at small 
q in the system modeled by Eq. ([!]). Eq. ([3| is the 
most important analytical result of this work and pro- 
vides a simple expression for the sound velocity that can 
be compared with on-going and future experiments on 
plasmons in G/M. In Fig. [I] we present a comparison be- 
tween the dispersion u ac (q) of the acoustic plasmon as 
found from the numerical solution of s(q,uj) = and the 
long- wavelength analytical result oj ac (q 0) = c s g, with 
c s given by Eq. (|3| (thin solid lines in Fig.[T]): the agree- 
ment is clearly excellent. 

The existence of acoustic plasmons has a large impact 
on the wave-vector and energy dependence of the one- 
body spectral function A(k, uS) [10] of G/M. As explained 
in Refs. [5] and [7] resonant interactions between charge 




FIG. 1: (Color online) Acoustic plasmon dispersion (in units 
of sf) in G/M. Filled circles and squares refer to different 
values of the graphene- metal gate distance d. Filled triangles 
refer to the standard plasmon (oc ^Jq for q —> 0) in the absence 
of a metal gate, here plotted only up to q = 0.4 /cf for clarity. 
These results are for \n\ = 5.0 x 10 13 cm -2 and a ee — 2.2. 
The long-dashed line (lj = vq) represents the upper bound 
of the intra-band electron-hole continuum. The intersection 
between the acoustic plasmon and the thin short-dashed line 
(lj = 2sf — vq, which is the lower bound of the inter-band 
electron-hole continuum) gives the critical wave number q c at 
which Landau damping starts. Note that, for every fixed q, 
the energy of the acoustic plasmon decreases with decreasing 
d, since the screening power of the metal becomes larger at 
smaller values of d. At large q the role of the metal gate is 
negligible and all plasmon curves collapse onto each other. 



carriers and plasmons give rise to strong plasmaron satel- 
lite bands, which are clearly visible in A(k,uS). In the 
system modeled by Eq. 0, resonances occur when the 
bare quasiparticle velocity v equals the acoustic-plasmon 
group velocity, i.e. when v — duj ac (q)/dq\ q=qi , = 0, at 
some specific wave number q*. If q* is smaller than 
the critical wave number q c at which inter-band Landau 
damping starts (see Fig. [I]), then a composite carrier- 
acoustic plasmon excitation, which we term for brevity 
"soundaron", emerges. Despite its simplicity, the model 
in Eq. ([!]) admits a non-trivial "phase diagram" for the 
existence of soundarons, which is shown in Fig. [2] We 
clearly see that soundarons exist in a very broad region of 
parameter space. When resonances occur, strong satel- 
lite bands appear in the one-particle Green's function 
A(k,uj). 

Quasiparticle decay rate and the one-body spectral 
function. — The spectral function A = A+ + A- 
can be written in terms of the quasiparticle self-energy 
£A(k,w) HQ]: 



A> 
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where £\(k) = Xvk — Sf are Dirac band energies. In 
the so-called a G W-RPA" approximation [5] [7J [8] [TO] the 
imaginary part of Y,\(k,uj) is given by the following ex- 
pression: 

9m [Z x (k,Lo)} = E / V *ti> ^~\QM 

x T x ,y{k,k+) [e(n) -e(-6'(fe + ))] 

(5) 

where k + = k + q, £1 = uj — <^/(fc + ), and J-^a' = 
[1 + AA / cos(6> fc)fc +)]/2. In Eq. ([5| 6(x) is the Heaviside 
step function, ^sm [£a] measures the band-quasiparticle 
decay rate. The two factors in the second line of Eq. ([5| 
express the influence of chirality and Fermi statistics on 
the decay process, respectively. The "quality" of the 
GoW-KPA approximation for doped graphene sheets has 
been carefully tested in Ref. [5] against angle-resolved pho- 
toemission spectroscopy (ARPES) data for graphene on 
SiC. Theory compares very well with the latter in a wide 
range of parameters [5], thus demonstrating that GqW- 
RPA is a very good starting point for doped graphene. 

Representative numerical results for [E\(k, uj)] and 
5fte [E^fcjCj)] in the presence of a grounded metal gate 
are collected in Fig. [3J The main features we want to 
highlight in these plots are the sharp resonances in the 
quasiparticle decay rates stemming from resonant carrier- 
acoustic plasmon interactions, as discussed above. Since 
5fte [E\(k,uj)] is related to [E\(k,uj)] by a Kramers- 
Kronig transform, sharp resonances in ^sm [E\(k,uj)] im- 
ply rapid changes in the real part of the quasiparticle self- 
energy and thus multiple solutions of the Dyson equation, 
uj — £fc,A — pA(fej^)] = 0- Those solutions to which 
a small ^sm [Ea(&,co>)] is associated give rise to sharp 
soundaron satellites in A(k, uj). In Fig. [4^) we present the 
dependence of the spectral function A(k,u) on momen- 
tum and energy (below the Fermi energy) for n-doped 
G/M. We clearly see soundaron satellite bands, which 
are especially strong close to the Dirac point (fe = 0). 
The dependence of A(k,ui) on a ee and d is illustrated in 
the supplementary material. For the sake of comparison, 
in Fig. [4|d) we have reported A(k, uj) for a graphene sheet 
with the same doping and at the same coupling constant 
(a ee = 2.2) but in the absence of a metal gate. Contrary 
to expectation, the plasmaron satellite bands in Fig. [4]o) 
are less pronounced than the soundaron satellite bands 
in G/M. This is in agreement with Ref. [5] in which it was 
clearly shown that plasmaron satellite bands become less 
intense and sharp at strong coupling (i.e. at large values 
of a ee ). 

In summary, we have shown that screened e-e interac- 
tions in a graphene sheet deposited on a metal gate are 
responsible for intriguing many-body effects. Despite the 
strong screening exerted by the metal, the electron gas 
in the graphene sheet displays self-sustained, long-lived 
oscillations which have a linear dispersion at long wave- 
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FIG. 2: (Color online) The color red (blue) labels the region of 
parameter space in which resonant interactions between car- 
riers and acoustic plasmons exist (do not exist) giving rise to 
"soundaron" satellite bands in the spectral function A(k,u). 
This phase diagram is for a coupling constant a ee — 2.2. 
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FIG. 3: (Color online) Quasiparticle decay rates and ener- 
gies for n-doped G/M. Panel a) |£sra[£+(fe, uj)] | (thick solid 
line) and Re [E+(fe, a;)] (dashed line) as functions of en- 
ergy uj (in units of and measured from the Fermi energy 
£f) for k = 0.1 far- The parameters used in these plots are 
n = 5.0 x 10 13 cm" 2 , d = 5 A, and a ee = 2.2. The in- 
tersections between uj — (dotted line) and Re [£+(&, a;)] 
(dashed line) represent the solutions of the Dyson equation, 
uj — — Re [£+(fc, uj)] = 0. Panel b) Same as in panel a) 
but for \%m[E-(k,uj)]\ and Re [£_(fc,u;)]. 



lengths. These modes yield composite carrier-acoustic 
plasmon quasiparticle satellite bands in the one-particle 
Green's function. Acoustic plasmons can be observed 
e.g. by inelastic light scattering or electron energy loss 
spectroscopy [14] [19]. Soundarons can be detected in 
tunneling density-of-states [20] and angle-resolved pho- 
toemission spectroscopy [5] [2T] experiments. 

While this manuscript was being finalized, we learned 
of an ab-initio study [22 that emphasizes the large im- 
pact of screening on the (high-energy) 7r plasmon of 
graphene. Our work is complementary to this study since 
we focus entirely on the low-energy "sheet" plasmon and, 
most importantly, on the spectral function. 
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FIG. 4: (Color online) Panel a) Spectral function A(k,uj) of 
n-doped G/M. Energy uj is in units of and measured from ep- 
This result is for n = 5.0x 10 13 cm -2 , d = 5 A, and a ee = 2.2. 
The solid lines represent the noninteracting bands. Panel b) 
Same as in panel a) but in the absence of the metal gate. The 
data in this panel have been obtained with the same theory 
and codes used in Refs. [5] and The only difference is that 
the results in these earlier articles are for smaller values of the 
coupling constant a ee (< 0.75). 
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SUPPLEMENTARY INFORMATION 

In this Section we report a series of numerical results 
to illustrate the dependence of the one-particle spectral 
function A(k,uj) on the graphene-metal gate distance d 
and fine- structure coupling constant a ee . For the sake of 



FIG. 5: (Color online) Same as in Fig. [4^l) of the main text 
but for d = 3 A. 



completeness, we also report data for the spectral func- 
tion of an isolated graphene sheet for different values of 

Qf ee . 

In Fig. [5] we show the GoW-RPA spectral function 
A(k,uo) of the electron gas in G/M for n = 5.0 x 
10 13 cm -2 , d = 3 A, and for a ee = 2.2. This result 
should be compared with that reported in Fig. [4^) of the 
main text, which was obtained by using a larger value 
of the graphene-metal substrate distance (i.e. d = 5 A). 
We clearly see that the soundaron satellite bands are ro- 
bust against changes in d. This is in agreement with the 
phase diagram presented in Fig. [2] of the main text: for 
large values of the carrier concentration (n > 10 13 cm -2 ) 
the critical distance d at which soundarons disappear is 
truly minute. 

In Fig. [6] we report plots of the GqW-KVA spectral 
function A(k,uS) of the electron gas in G/M for n = 
5.0 x 10 13 cm -2 , d = 3 A, and two values of the coupling 
constant a ee . Comparing these results with those re- 
ported in Fig. [5j we clearly see that the soundaron satel- 
lites are also robust against changes in the coupling con- 
stant. Notice, in particular, that intensity and sharpness 
of the soundaron peaks increase with decreasing a ee - Of 
course, in the noninteracting a ee — >• limit these bands 
smoothly disappear. 

In Fig. [7] we illustrate the dependence of the spectral 
function of an isolated (i.e. in the absence of a metal 
gate) graphene sheet on a ee . These results should be 
compared with those reported in Fig. 4b) of the main 
text. Also in this case, the satellite bands due to resonant 
carrier-plasmon interactions are sharp at moderate-to- 
intermediate values of a ee . This is in agreement with 
the results reported in the supporting online material of 
Ref.EJ 
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FIG. 6: (Color online) Same as in Fig. [5] but for two smaller FIG. 7: (Color online) Same as in Fig.JlJ)) but for two smaller 
values of the coupling constant a ee - values of the coupling constant a ee - 
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